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1. Introduction 

Quantum Chromodynamics (QCD) is the theory of interactions between quarks and gluons. 
As is well known, this theory is asymptotically free at high energies, while at low energies, con- 
finement sets in. There is still no good understanding of confinement as non-perturbative aspects 
start to play an important role at the relevant energy scale. Therefore, it is of high interest to study 
pure gluonic particles where only gauge fields play a role. These pure gluonic particles are called 
glueballs. 

So far, there has been 20 years of intensive experimental search towards glueballs, see ^ for a 
complete overview on the experimental status. Unfortunately, no definite answer to the question 
whether a glueball has been observed or not can be given. Although, in the scalar sector, there is 
no question that more states have been found than can be accommodated by a single meson nonet. 
Hopefully, the future will provide us will more data on this as there are many new experiments 
planned, e.g. the PANDA Experiment at GSI in Germany [§], BES III at BEPCII in Bejing [§], the 
GlueX Experiment at Jefferson Laboratory in the USA [§], ALICE at CERN [§]. 
As no clear experimental data is available on glueballs, lattice QCD has become very important for 
extracting relevant information about glueballs. In the quenched approximation, the state of the art 
is currently given by [||, [7|]. For example, the masses of the three lowest lying glueballs are given 
by 

scalar glueball : ++ = 1.710GeV/c 2 , 
tensor glueball : 2 ++ = 2.390GeV/c 2 , 
pseudo scalar : + ~ = 2.560GeV/c 2 . (1.1) 

On the theoretical side, various models and methods have been developed to calculate properties 
and masses of the glueballs. Historically, the first model to study glueballs is called the MIT bag 
model, and was applied to glueballs in 1976, see After this, many other models have been 
constructed, often phenomenological or with at least some input from experimental data. For an 
extensive overview, see [§]. 

A very fundamental way to investigate bound states of gluons, is to study the correlator 
(Fjiv( x )Fap(y))' where F 2 V is classically the gauge field strength. In the next sections, we shall 
study this correlator using first the Faddeev-Popov action, and then the much more complicated 
Gribov-Zwanziger framework. In both theories, we shall show that one can construct a renor- 
malization group invariant containing F 2 V . However, the Yang-Mills gauge theory obeys a BRST 
symmetry s, while the Gribov-Zwanziger action breaks this BRST symmetry. This has big conse- 
quences for the correlator, as we shall show in the next sections. 

2. The construction of the glueball operator in the Yang-Mills action 

2.1 The Yang-Mills action and the BRST symmetry 

When turning to the quantumlevel, one can overcome the problem of integrating over zero 
modes in the Yang -Mills action J d 4 x|F 2 v , by introducing a gauge fixing. When choosing the 
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Landau gauge, d^A a ^ = 0, the partition function is given by 

Z FP = J [dA]8(dA) teLM ab e-if* xF Z* F Z* , (2. 1) 

with 

M ab = -d^d ab -gf bc A c li ), (2.2) 

the Faddeev-Popov operator. We can lift the determinant into the action, by introducing the ghost 
fields c and c, 

Z FP = J [dA] [dc] [dc] [db}e~ SFP , (2.3) 

with 

Sfp = \J d 4 xF« v F° v + J d 4 x (V<V« + c a d^D ab c b ) . (2.4) 
We have also introduced the auxiliary ft-field, which implements the Landau gauge. 

The Faddeev-Popov action is invariant under the following nilpotent BRST transformation, 



with 



sS FP = , s 2 = , (2.5) 



sA^ = - (D^c) a , sc a = \gf abc c b c c , sc a = b a , sb a = . (2.6) 



2 

This BRST symmetry is of paramount importance as it is at the origin of the Slavnov-Taylor iden 



tity, which enables us to prove the renormalizability of the Yang-Mills action [10, 11]. 



Another important property is that the BRST charge allows us to to define the sub-space of the 
physical states: 

sO phys = , (2.7) 

whereby O p hy S + sX is in the same class of operators as O p h ys , since s 2 X = trivially. The BRST 
symmetry s also allows us to establish the unitarity the S matrix at the end [fi^, |13[|. 



2.2 The construction of a physical operator 



Now, we can construct a scalar glueball operator containing F 2 which is renormalizable rfTj]. 



For this, we introduce a dimensionless source q to couple F 2 V to the Faddeev-Popov action Spp. 
Naively, one would expect the following action to be renormalizable, 

S FP + j d A xq{x)F 2 v (x) . (2.8) 



However, F 2 V will mix with other dimension four operators and therefore the action (2J3) is not 
renormalizable. 

There exist three types of dimension four operators. Firstly, we have the gauge invariant opera- 
tors J? '. These are constructed from the field strength F^ v and the covariant derivative D a ^. These 
are BRST closed but not exact, e.g. F 2 V . Secondly, we have the gauge invariant exact operators S, 
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e.g. s^dfiA"^), with s the BRST variation (2.6). Finally, one can construct "Equation of motion 



terms" J>t, e.g. A° which vanish upon imposing the equation of motion. 

Now, one can intuitively easily understand that these 3 different classes will mix in a certain way 



[ ]15[ , |16| , |j7Q . Firstly, a <§ operator cannot mix with a gauge invariant operator & ' . Indeed, as 

(<f ) = (*(•••)) =fj, (2.9) 



and suppose $ mixes with the two other classes of operators, 

{S Q )=a{^)+b{S) +c{Jif) 
=o =o 



(2.10) 



we know that the mixing coefficient a has to be equal to zero. Secondly, a J4? operator cannot mix 
with & and $ as <ffl will vanish upon using the equation of motion, while $ as ffl will not. This 
can be translated in a upper triangular form of the mixing matrix, 




/ Z&g Z&g 
Zgg Zgjtf 

\ o z^^ 







) 





(2.11) 



The arguments for this mixing behavior are very general, and therefore expression ( |2. 1 1[ ) should 
be valid for every gauge theory which has a BRST symmetry s. 



With this knowledge, we can construct an improved action of ( |2.8| ) which is renormalizable, by 
introducing other dimension 4 operators. We shall couple a source to all of the following operators 



1 



4 w ' 
• = s^d^) , 
8Sfp 



J% = At, 



(2.12) 



and analyze the renormalization following the formalism of algebraic renormalization [JTCj] . Details 



of the calculations can be found in [14]. We have proven that the new action including the three 



operators in expression (2.12) is renormalizable, whereby the mixing matrix is given by 









ZjqZ qq 




\- 





1 






V 






-Zi„Z, 



Jq^qq 


1 



(2.13) 



Notice the upper triangular form predicted in expression (2.1 1 ), which is now clearly proven. We 
also see that $ does not mix with M 3 , which is logical as $ is a sum of two "equation of motion" 
terms and behaves like J^. 
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Now that we have found the mixing matrix, we shall completely fix the entries of this matrix 



following the lines of []18fl. Starting with the following most general (n + 2m + r) -point functions 

yn + 2m + r = (^[xi) . , .A(x n )c{y\) . ■ ■c(y m )c(y l )c(y m )b(zi) ■ • •£(■£/•)) 

[d&]A(xi). . .A(x n )c( yx ). . .c{y m )c{h)c{y m )b{z,i) . ..b{z r )^ S , 



(2.14) 



it is possible to determine the mixing matrix (2.13) to all orders by using the fact that dB 



be completely finite. We have found that [14] 



dg 1 
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1 / 




(2.15) 



whereby j3 is the /3 -function and y c is the anomalous dimension of the ghost. Let us mention that 
this latter technique of determining the mixing matrix to all orders is only possible for the operator 
mixing matrix of the operator F 2 V . This would not work for other operators e.g. F\F^ a F va — 
^8[i V F 2 , etc. This is related to the fact that deriving the Green's functions (2.14) w.r.t. g 2 brings 
down the operator of interest F 2 V = ^ from e~ s . 



The final step is to determine a renormalization group invariant from fl2.15| ). This can be done 
easily and we have found, 



phys 



(2.16) 



This means that instead of (^(x)^(y)}, we should actually investigate the quantum version 
(&{x)M(y)). We can make a few observations. Firstly, we see that this renormalization group 
invariant 2% coincides with the trace anomaly @{J , which is expected as 0^ is also a d = 4 renor- 



malization group invariant [19]. Also, we recover the same mixing as in some older works in a 



slightly different formalism, see [g(], 2TJ . Secondly, we can immediately set the term in equal 
to zero as in the correlator we are working on-shell. Finally, the BRST symmetry plays an impor- 
tant role in the correlator. Due to s-invariance of the action, we find 



(O phys (x)O phys (y)) = {&(x)&(y)) 



4*2 V +*(-..) 



to 



F£ v + s(. 



(v) 



|r) F 2 v (x)F 2 v (y) + s(...) 



[F^ v (x)F 2 v (y)) 



(2.17) 



meaning that the mixing can be ignored in the Yang-Mills action at the level of the correlator due 
to the BRST symmetry. 
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3. The construction of the glueball operator in the Gribov-Zwanziger action 

3.1 The Gribov-Zwanziger action 

In 1977, Gribov discovered that the Landau gauge is plagued by the existence of Gribov copies 



[22], i.e. there are still multiple gauge copies which all fulfill the Landau gauge condition. We 
can get partly rid of these copies by restricting the domain of integration in the Feynman path 
integral to the Gribov region Q. [p3|], 

D. = {A;dA = 0,-dD>0} , (3.1) 

or thus by demanding the Faddeev-Popov operator to be strictly positive. The partition function 
( ^T| ) is then translated into 

Z H = J [dA]8(dA) detM u V(^ d4v Wv+ 5 ") , (3.2) 

with Sh a non-local term 

S H = y A g 2 J f^A^ {M- l ) ad 'f dec ; A* , (3.3) 



the so-called horizon term Q23|]. As a non-locality is difficult to handle, it is useful to recall that one 



can cast this action into a local form 

S GZ = S F p + So + S r , (3.4) 

with 



So = Jd A x(lp a jM ab (p* c -ld a jM ah Q) b ^ 
S r = -fg f d 4 x (f^^+^K + - g {N 2 - 1) y 2 ) . (3.5) 



We have introduced a pair of complex conjugate bosonic fields y(p a ^, (p^ j and a pair of Grassmann 

fields (co a ^ , co^ c ). We recall that Sqz is renormalizable to all orders [24, 25, 26], even in the pres- 
ence of massless [27, 28] or massive quarks [pSQ. 



The parameter y is called the Gribov parameter, has dimension of a mass and is not free but deter- 
mined by 



dy 2 



^3=0, (3.6) 



with r the quantum action defined as e r = J[D<J>]e s . 

Let us have a closer look at the BRST breaking of the Gribov-Zwanziger action. The new fields 
have the following BRST transformation: 

j<pf = fi)f, sw-=0, sw?=-$, #-=0. (3.7) 
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One can check very easily that Sgz is no longer invariant under the BRST symmetry s, 



sS GZ = s{S YM + S gf + S» + S r )= S {S r )=gY 1 J d d xf bc farf (<p» c + <pj e ))(3.8) 



This soft breaking has been studied in Q33Q. Let us remark that despite this breaking, the Gribov- 
Zwanziger action is renormalizable, due to a rich set of Ward identities. Moreover, only two renor- 
malization constants are needed, Za and Z g , just as in the Faddeev-Popov action. As a byproduct of 
this action, we see that the gluons get "confined" by the horizon. Indeed, as the gluon propagator 
has complex poles, gluons cannot describe physical excitations. In addition, we also observe that 



the Fourier transform shows positivity violation [33]. The question which rises now, is how to 



define physical operators as the definition (2.7) is no longer possible here 



Before giving a possible answer to this question, let us first mention that it is possible to embed 
Sqz into a "larger" action. For this, we replace Sy with S s = s{. . .), 

Sgz = SyM + Sgf + So + Sy 

I 

£gz = Sym + Sgf + So + S s 

We have introduced new sources in S s , but in the end, we can set these sources to the right values 
so that S.s|pi,y S = Sy. This implies that we did not change the theory. The advantage is that we have 
an action which is BRST invariant again. 



3.2 The construction of a potentially physical operator M 

With £gz we can construct a glueball operator which is renormalizable, completely analogous 
as in the Yang-Mills case, and in the end we can take the physical limit to the physical Gribov- 



Zwanziger action Sqz- This is a highly non-trivial task, see [30] for all the details. However, the 



final outcome is very peculiar. We find exactly the same expression as in (2.16), 



Ophys — & 



' -2y c > 



Iphys ^Yc 



3%" 



phys ' 



(3.9) 



but $' is no longer BRST exact: 



Iphys 



3t 



s(d^ + d^fDtftf) + fDf (<p^ + <pj*j + d(N 2 - 1) 7 4 , 
<SS G z 



phys 



A?, 



(3.10) 



The fact that expression ( |3.9| ) has exactly the same form as ( |2.16| ) is closely related to the fact that 
the Gribov-Zwanziger action has only two independent renormalization constants, Za and Z g just 
in the ordinary Yang-Mills action. 



This time, $' will not be trivial at the level of the correlator. Not only is $' no longer BRST 
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invariant, also the action Sgz breaks BRST softly. Therefore, 

(O phys (x)O phys (y)) = im\x)M'iy)) 



(x) 



r^-2y c <f 



(v) 



(3.11) 



and thus the correlator has extra parts in comparison with the standard glueball correlator 

j8 x 2 



JF(x)JP(y)). In fact, Zwanziger [23] has studies the pole structure of 



d 4 x(F 2 v (x)F 2 v (0)) e ^ 



(3.12) 



at tree level and he has reported that this correlator displays a physical cut at p 2 = —ly 2 , while it 
also displays unphysical cuts at p 2 = ±4iy 2 . It would therefore be interesting to investigate the cut 
structure of the improved correlator (M' (x)& ! (y)) , to see whether it can be called physical or not. 

3.3 An extra word about contact terms 



There is perhaps one point which needs to be clarified w.r.t. contact terms. In [ 14, 3C] we have 



always safely neglected contact terms as such terms are usually ignored Q15Q. The origin of such 
contact terms is the following. The correlator (f 2 v (x)F 2 p (y) J stems from deriving the partition 
function w.r.t. q(x) and q(y) and then setting all sources equal to zero. One can imagine, having a 
pure source term like / d 4 q(x)d 4 q(x) into the action, this gives rise to a contact term ~ d 4 8(x — y). 
However, as we are not interested in the space time point x = y, this term shall always be zero, since 
5(0) =0 in dimensional regularization. Therefore, we have always neglected this class of terms. 



One can object now, as going to Fourierspace, this term is no longer irrelevant. Indeed, 

d 4 8{x-y) = J-jtotffM (3.13) 

giving rise to a term ~ k 4 . It is therefore possible that a term like this will be necessary to kill 
divergences of the form y, present in the Fourier transform of (^F 2 v (x)F 2 p(y)^ . As already said, 
when going back to x- space, a term like this is zero in dimensional regularization for x^y. More- 
over, in Fourierspace, contact terms will only generate polynomial contributions in the momentum 
and/or external mass scales. Such terms are therefore irrelevant when investigating the pole and cut 
structure of the propagator in momentum space. 



In conclusion, contact terms are irrelevant for the calculation of the pole and cut structure of 
(M' (x)M' (y)) , as they do not contribute. They are however needed to ensure a completely finite 
correlator function, see [Il5|, [311]. 



4. Conclusion 

We have investigated the renormalization of F 2 V using the standard Faddeev-Popov action, as 
well as the more complicated Gribov-Zwanziger action, which takes into account Gribov copies. 
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Firstly, in the Yang-Mills case, we have found that F 2 V mixes with a BRST exact operator S 
s(c a d^A a ^). We have constructed a renormalization group invariant 



!^l^-2y c (g 2 )<? . (4.1) 



However, at the level of the correlator, this BRST exact operator does not play any role, 

(M(x)M(y)) = (^p) {Fl v {x)F^ v {y)) , (4.2) 
due to the s-invariance of the Faddeev-Popov action. 

Secondly, in the Gribov-Zwanziger case, something peculiar happens. The issue of the Gribov 
copies leads to a modification of the Faddeev-Popov formula. The BRST invariance of this modi- 
fied Faddeev-Popov action turns out to be softly broken. After studying the renormalization of F 2 V , 
we have found again that F 2 V mixes with <?' = s(c a d^ + dffiDfiqt*) + fDf(<p b J + 9*?) + 
d{N 2 -\)f. Notice that g' is no longer BRST exact. The renormalization group invariant has 
the same form as in the Yang-Mills case, 

M' = ^-&-2y c {g 2 )g' . (4.3) 
r 

However, now g' does not vanish at the level of the correlator, as the BRST is softly broken in the 
Gribov-Zwanziger action. 

{M{x)M(y)) = {^r) (&{xW(y))-2 yc ^(^(x)g(y)) 

-2 7c 4 (g(x)F 2 (y)) + 4 7 2 (g(x)g(y)) . (4.4) 

The possibility that the breaking of the BRST in the Gribov-Zwanziger action might be relevant 
in order to construct the local physical operators in the presence of the Gribov horizon looks very 
attractive. This will be investigated in future work. 



Finally, as a last remark, let us mention that this analysis is not affected by taking into account 
dimension two condensates into the theory. Therefore, we would find exactly the same results 
using the so-called Refined Gribov-Zwanziger framework, see p2 , [33]] . 
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